Pseudogaps in the t-J model: Extended DMFT study 
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We investigate the highly incoherent regime of hole-doped 2d Mott-Hubbard insulators at mod- 
erately small doping S and temperatures > O.IJ, where J is the exchange coupling. Within an ex- 
tended dynamical mean-field theory of the t-J model and a generalized non-crossing approximation 
we calculate the single-particle spectral function, the dynamical susceptibility, and thermodynamic 
and transport quantities. Short-ranged antiferromagnetic fluctuations lead to strongly incoherent 
single-particle dynamics, large entropy and large electrical resistivity. At low doping a pseudogap 
is found to open up both in the single-particle and the spin excitation spectra leading to a decrease 
in entropy and resistivity. The Ifall coefficient changes sign to positive values upon lowering the 
doping level and increases inversely proportional to S. 

PACS numbers: 71.30.-|-h,74.72.-h,71.10.Hf 
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I. INTRODUCTION 

Strongly interacting Fermi systems on a two- 
dimensional lattice have been a focus of interest ever since 
high temperature superconductor materials were discov- 
ered in 1986. The unusual properties of these materials, 
in particular in the normal conducting phase, have led 
to the hypothesis that the usual Landau Fermi liquid 
theory may not be applicable in this casei. A widely 
accepted view holds that these systems may be consid- 
ered as hole-doped Mott-Hubbard insulators. The corre- 
lations present in a Mott-Hubbard insulator are charac- 
terized by strongly suppressed charge fluctuations (due 
to the constrained hopping of the holes resulting from 
the strong on-site Coulomb repulsion U) and enhanced 
quantum spin fluctuations governed by the antiferromag- 
netic nearest-neighbor spin exchange interaction. The 
interplay of the motion of holes with the antiferromag- 
netically correlated spin background is the central prob- 
lem of hole-doped Mott-Hubbard insulators. Despite an 
extraordinary effort by many theorists and a correspond- 
ingly large number of papers we believe it to be fair to 
say that a thorough understanding of this problem is still 
lacking. 

The ground state of the Mott-insulating state of elec- 
trons on a square lattice at half-filling is expected to 
be antiferromagnetically ordered^. Doping with holes 
leads to a rapid destruction of long-range order, at a 
critical concentration Sc of a few percent doping. For 
larger dopings there is evidence for strong antiferromag- 
netic spin fluctuations of relatively short range. 

In this paper we undertake to explore the consequences 
of strong incoherent and local (i.e. nearest neighbor) spin 
fluctuations on the dynamics of charge carriers, and on 



the thermodynamics of the systems. We also investi- 
gate how the single-particle properties feed back into the 
spin dynamics. Our approach is focused on the temper- 
ature regime of O.IJ ^ T < t {J = 0.3t in the cuprates) 
where J is the exchange constant and t is the nearest 
neighbor hopping amplitude. In this regime we expect 
strong quantum and thermal fluctuations driven by com- 
peting interactions to decohere the fermionic excitations. 
This temperature regime is bounded from below by possi- 
ble antiferromagnetic, superconducting or other ordered 
states. The incoherent regime is confined to small doping 
levels Sc ^ S ^ 0.3, and crosses over into a Fermi liquid 
state at 6 > 0.3. 

A minimal model encompassing the physics described 
above is the t-J model. It is well known that the in- 
terplay of hopping and the local correlations induced by 
the on-site Coulomb interaction may be captured in Dy- 
namical Mean-Field Theory, in which the lattice model is 
mapped onto a quantum impurity coupled to a fermionic 
bath in a self-consistent fashion^. In the same spirit 
the nearest neighbor exchange interaction of a given spin 
to its neighboring spins may be approximated by a dy- 
namically fiuctuating bosonic field, to be determined self- 
consistently^-^. In this way the two principal processes, 
constrained hopping and spin exchange interaction, may 
be fully incorporated on the same footing, on the level of 
short-range correlations. 

As reviewed in Section |h1 the Extended Dynamical 
Mean Field Theory (EDMFT) for the two-dimensional t- 
J model is obtained by approximating the single-particle 
self-energy Sk(t^) and the two-particle self-energy Mq(w) 
by momentum independent functions. S(ct;) and M{ll!) 
are obtained by equating the local (i.e. the momentum 
integrated) single-particle Green's function and spin sus- 
ceptibility, respectively, with the corresponding quanti- 



2 



ties of an extended Anderson impurity model, featuring 
a fermionic and a bosonic bath to be determined self- 
consistently. For the bare hopping integrals and exchange 
couplings we use a nearest-neighbor tight-binding model 
on the square lattice. The local approximation is bet- 
ter the higher the spatial dimension d and becomes ex- 
act for d — > oo, provided the hopping amplitude t and 
the exchange coupling J are scaled as t/^/d and J/Vd. 
This scaling is possible in the paramagnetic regime. Most 
of the methods employed for the solution of the Ander- 
son impurity or Kondo problem do not work here. We 
use self-consistent perturbation theory in the form of 
conserving approximations^-, and the exact projection 
onto the Hilbert space without double occupancy (limit 
U oo). We are interested in describing the highly in- 
coherent regime at small doping levels and not too low 
temperatures, where the spectral functions are broad and 
relatively featureless. In this regime we expect vertex cor- 
rections and higher order processes in general, to change 
the characteristic parameters like maximum values, peak 
widths and gap widths of the dynamic quantities by cor- 
rection terms of order unity, but we do not expect that 
these contributions lead to more coherence or new col- 
lective behavior. In this spirit we approximate all self- 
energies by their lowest order self-consistent perturba- 
tion theory expressions (in the hopping parameter and 
exchange coupling). The resulting theory, presented in 
Section imi is an extension of the Non-Crossing Approx- 
imation (NCA)1S, including the bosonic bath. 

The results of this approximation scheme for the t-J 
model are presented in Section IIVI It turns out that 
nearest-neighbor spin fluctuations are sufficient to cre- 
ate a pseudogap in the single-particle spectrum and in 
the spin excitation spectrum at g-vectors away from 
(tt, tt), for small dopings 6 < 0.1, similar to what 
is seen in ARPES experimentsii and in the magnetic 
propertiesiSii^. The pseudogap scales with J. There are 
several indications that Fermi liquid behavior is violated 
for S < 0.2. Most noteworthy, the effective chemical po- 
tential is found to move from the center of the band up 
to the band edge, as the doping is decreased to small 
values. As S grows beyond 0.25, however, Fermi liquid 
behavior appears to be restored. The entropy turns out 
to be large in the range 0.1 < S < 0.2 and is reduced 
on both sides of this interval by the pseudogap and the 
incipient Fermi liquid behavior, respectively. The resis- 
tivity is dominated by strong incoherent scattering, and 
the Hall coefficient is found to be hole-like, c>c 1/d, for 
small S, again resembling the observed behavioi»i4. Some 
of the results have been reported in • 

Results similar to ours have been found in two re- 
cent works using DMFT for a clusters of sites within 
the Hubbard model. Maier et alm^ applied the dynami- 
cal cluster approximation (DCA) for various cluster sizes 
up to 64 sites to the Hubbard model in the intermedi- 
ate coupling regime {U ~ bandwidth). The DCA equa- 
tions were solved with Quantum Monte Carlo techni ques 
down to room temperature. The authors of Ref. [ig 



identified signals for non-Fermi liquid behavior at low 
doping i5 < 0.1 and found a large residual scattering rate 
and a pronouced pseudogap at low doping. In Ref. [T^ . 
Stanescu and Phillips studied the Hubbard model at in- 
termediate coupling within a two-site cluster approach 
using the Non-Crossing Approximation as a quantum im- 
purity solver at not too low temperatures. It is again 
found that Luttinger's theorem appears to be violated 
for low doping in a regime where a pseudogap opens. 

Despite the similarity of the numerical results, quite 
different explanations for the observed pseudogap physics 
have been suggested, ranging from short-range spin cor- 
relations, spin-charge separation and resonance valence 
bond (RVB) physicsi^ to effects of the upper Hubbard 
band and current-correlations involving three neighbor- 
ing sites^^. By construction our approximation scheme 
is not able to describe such intersite correlations or 
RVB singlets and does not include the upper Hubbard 
band: nevertheless the overall results are qualitatively 
very similar. We take this as a strong indication that 
neither short-range magnetic or current correlations nor 
RVB physics is the underlying reason but argue that 
there is another generic mechanism for pseudogap for- 
mation: The strongly incoherent dynamics captured in 
our scheme as well as those of Refs. [ma appears to 
be the dominant feature of the Hubbard model as well as 
the t-J model in the low energy sector (O.IJ ^ w < i) for 
small doping. Therefore pseudogap formation seems to 
be a generic property of any strongly incoherent Fermi 
system close to a Mott insulator. In other words, the 
existence of a pseudogap neither requires slowly fluctu- 
ating, finite-ranged ordered domains (antiferromagnetic, 
superconducting)^^ nor a local resonance state. 



II. EXTENDED DYNAMICAL MEAN-FIELD 
THEORY OF THE t-J MODEL 

The standard model embodying the physics of the hole- 
doped Mott-Hubbard insulator is the t-J model, defined 
by the Hamiltonian 

where Si = ^ X^o- <t' ^^''o-o-'Cicr' is the spin operator at 
lattice site i, r denotes the vector of Pauli matrices 
and tij(Jij) are the hopping amplitudes (exchange in- 
teraction) connecting sites i and j. For the numeri- 
cal evaluation to be discussed later we will use a tight- 
binding model on a two-dimensional square lattice, tij = 
—tSi^i+r, Jij — JSi^i+T, where r labels nearest neighbor 
sites. The operator c^yicia) creates (annihilates) an elec- 
tron at site i with spin projection cr at a singly occupied 
lattice site. In terms of usual electron operators cf^{cicr) 
one has c+ = c+ (1 - n^,-a), where riia- = cf^Ci^ is the 
occupation number operator. In this way occupation of 
lattice sites by two electrons with spins "f and | is avoided. 
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which would cost the large Hubbard energy U . We will 
be interested in electron densities close to half-filling of 
the band, such that (X^o- '^ic) — n — 1 ~ 5, where (5 <C f 
is the doping concentration of holes. 

Whereas at exactly half-filling, when H reduces to the 
Heisenberg model, the ground state has antiferromag- 
netic long-range order, we anticipate that this will not 
be the case for sufficiently large doping 5 > 5c (in ex- 
periment 5c — 0.03). In this regime it is reasonable to 
assume the antiferromagnetic correlations in the system 
to be short-ranged. We assume furthermore that addi- 
tional forms of long-range order (such as superconductiv- 
ity) that may be possible ground states of the t-J model 
are confined to a lower temperature regime, such that 
the corresponding fluctuations are sub-dominant at el- 
evated temperatures. Consequently, one expects an ex- 
tended high-temperature regime where short-ranged spin 
fluctuations lead to a highly incoherent metallic state, 
as observed in high-temperature superconductors, with 
anomalous transport properties (large, non-Fermi liquid 
type electrical resistivity, hole- like Hall constant), large 
entropy, broad "quasiparticle" peaks in photoemission, 
etc. It is our aim to investigate this regime within an 
approximation scheme which neglects most of the longer 
range spatial correlations, but keeps the dominant short 
range spin correlations. 

The single-particle dynamics and the two-particle dy- 
namics of the model are described by the Green's function 

Gk,.(ic^) = - / rfre'""(?^r5k.(r)5ka(0)) 
Jo 



io; -I- /X — Ek — Sk,(T(iw) 
and by the spin susceptibility 

X^,c.m - / dTe*''M7^r5_q,„(T)5q,„(0)) 
Jo 



Here /3 is the inverse temperature T (we employ units 
with ks = h = I), uj and f2 are fermionic and bosonic 
Matsubara frequencies, and ek and Jq are the lattice 
Fourier transforms of the hopping amplitudes tij and 
the exchange couplings Jy, respectively. While the self- 
energies Yi\^{iLo) and Mq(zu;) are momentum dependent 
in general, the observation that the fluctuations in the 
system are short-ranged in the regime we are interested 
in suggests that a "local" approximation, neglecting the 
momentum dependence of E and M altogether, may be 
a good starting point. We therefore employ in this paper 
the main approximation 

Ek(w) ^ S(a;) (4) 

and 



thus capturing the effect of local fluctuations in time, 
which we expect to be important in the presence of strong 
inelastic scattering. 

The momentum independence of S and M allows us to 
map the lattice problem onto an Anderson impurity prob- 
lem where the host medium has to be determined self- 
consistently. Considering first the single-particle prop- 
erties, i.e. S(a;), the corresponding Dynamic Mean- 
Field Theory (DMFT) has been widely used to calculate 
properties of the Hubbard model and periodic Anderson 
model44&. One maps the problem onto an Anderson im- 
purity embedded in a fermionic bath. Applied to the t-J 
model it amounts to treating the exchange interaction in 
mean-field theory. This is not sufficient to allow us to 
maintain the balance between dynamical hopping pro- 
cesses and spin fluctuations, which is at the heart of the 
t-J model. We therefore follow Refs. 03 and extend the 
dynamical mean-field idea for the paramagnetic phase by 
introducing a fluctuating magnetic field coupling to the 
local spin as representing an additional class of degrees 
of freedom of the medium. This type of approximation, 
termed "extended DMFT" (EDMFT), has been appHed 
to the Kondo lattice model^ and the extended Hubbard 
modeli. It is important to note that EDMFT becomes 
exact in the limit of infinite dimensions d — > cx) , provided 
t and J are scaled as t/^/d and J/^/d, respectively. We 
shall use this property in deriving the EDMFT equations 
(see Appendix . We will, however, regard EDMFT as 
an approximation applied in finite dimensions, and as 
such will use the tight-binding expressions for ek and Jq 
valid in d = 2. 

To summarize, the EDMFT is probably best visualized 
by considering a single-site, the "impurity" , and its cou- 
pling to the surrounding "medium" . There are two types 
of coupling processes, as is evident from the Hamiltonian: 

(i) hopping to and from the "impurity" into the medium, 
as in the Anderson impurity model (in the limit of infi- 
nite t/, as a consequence of the no double occupancy 
constraint) . The medium is modeled by a non- interacting 
fermion system (the "conduction electrons" ) , whose local 
density of states has to be determined self-consistently. 

(ii) exchange coupling of the local spin at the "impurity" 
site to the spins of the medium. In the limit d — > oo 
the two components of the medium, fermions (see above) 
and spin fluctuations are completely decoupled. We do 
not expect that this approximation holds in 2d for low 
temperatures. But in the regime considered in this pa- 
per, where electrons are highly incoherent, we believe 
that such a modeling is appropriate. The spin fluctua- 
tions of the medium is described by a (vector) bosonic 
bath, whose spectrum again has to be determined self- 
consistently. 



Mq(w) ~ M{lo) , 



(5) 



In this way one is led to a generalized quantum impurity 
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model with Hamiltonian: 

ffEDMFT = ^Ekcl^Cka +V'^{cl^da + h.C.) 

k(T ka 

-find + Y.'^M • + Srf • (h, + hi,). (6) 

A formal derivation of i?EDMFT in the limit c? — s- cx) is 
given in Appendix ^ Here djr is a projected fermion 
creation operator for the impurity orbital (the original 
operator c^^ at the chosen "impurity" site 0), Ud — 
J^a'^tda and = ^J2a,a' '^t'^<^<^'d<^' ■ ^hc fermionic 
bath is represented by free fermion operators cj^^, the 
bosonic bath by free boson operators /ij^, a = 1,2,3 

with = {hqi,hg2,hg3), and J^qO^g + ^-q) playing the 
role of a fluctuating local magnetic field. The excita- 
tion spectrum of the bath degrees of freedom, Ek and 
Uq, as well as the coupling constants V and / have to be 
determined self-consistently by equating both the single- 
particle Green's function and the spin susceptibility of 
the impurity model Gimp,Ximp with the local Green's 
function Gioc and the local susceptibility xioc of the lat- 
tice model. 



Jo 

The local G and x are obtained from their lattice coun- 
terparts ®, taking into account Q), |(SJ), and by 
summation over all momenta 



Gioc{uj) = ^ Gk{iuj) 

k 



(8) 



(9) 



As shown in Appendix ^ the self-energies S and M 
also characterize the impurity Green's functions: 

Gimp(«w) = [io; + ^ - V^^Gc(ic^) - S(ic^)] ^ , (10) 



where 



X^n.p{ii0) = [M^I^Gh] 



iuj — Ek ' 



2u>„ 



(11) 



(12) 



(13) 



so that the system of equations |5J - (lll|l is closed. It 
follows from (|10I 10') that only the densities of states of 
the baths, 

Ac{uj) = — Im Gc(c^ - tO) ^V^y S{uj - Ek) (14) 



and 



Dh{uj) = — Im Gh{^-iQ) 

T 



E ['5(^ - "^q) - + ^9) 

1 



(15) 



are needed. For practical purposes we have included the 
coupling constants V and /, respectively, in the defini- 
tions of the density of states. 



III. GENERALIZED NON-CROSSING 
APPROXIMATION 

The solution of the quantum impurity model for 
given Aciio) and I?/i(aj) is difficult. Many of the meth- 
ods developed in the past for solving Anderson impurity 
models in the context of DMFT such as iterated pertur- 
bation theory^ and the numerical renormalization group 
methodic are not applicable in the case of a bosonic bath. 
The Quantum Monte Carlo method has been successfully 
applied to an anisotropic Kondo lattice model with Ising- 
type spin couplingSSiSi, but it is extremely difficult to 
treat Heisenberg couplings with manageable effort. The 
only method left to us is self-consistent perturbation the- 
ory like the non-crossing approximation (NCA) or the 
conserving T-matrix approximation (CTMA)24i&. 

We will therefore employ a conserving diagrammatic 
approximation in which infinite classes of perturbation 
theory in V and / are resummed. We are aiming at a level 
of approximation corresponding to NCA for the usual 
Anderson model. A convenient way to phrase the per- 
turbation theory in the hopping V and the exchange cou- 
pling J, in the presence of an infinitely strong Coulomb 
repulsion ?7, is in terms of a pseudo-particle representa- 
tion. We define pseudo-fermion operators f^,(7 =T,i, 
creating the singly occupied impurity state and the slave 
boson operator 6"*" creating the empty impurity level, 
when acting on a corresponding vacuum state~. Since 
the local level is either empty or singly occupied, the op- 
erator constraint Q = b^b + J2a fal<y = 1 ^as to be 
satisfied at all times. The constraint is enforced exactly 
by adding a term XQ to the Hamiltonian and taking the 
limit A ^ 00 (see (I2r>|l below). The projected local elec- 
tron operators d^ may then be replaced by b^ fa, turning 
the problem into a many-body system of pseudo-fermions 
fa- and slave bosons 5, interacting with the fermions Cka 
and bosons of the bath. 

It is essential for any approximation scheme to re- 
spect the projection and not to allow transitions be- 
tween different sectors of Hilbert space labeled by Q. To 
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a) 



b) 



<I> = 



FIG. 1: The two lowest order contributions to the Luttinger- 
Ward functional $ and corresponding self-energies. Only di- 
agrams with no line- crossings are taken into account (a gen- 
eralization of NCA). The broken (wavy) line denotes pseudo 
fermion (pseudo-boson) Green's function G/ (Gb), and the 
solid lines represent the conduction electron Green's functions 
Gc, the curly line the correlator Gh of the bosonic bath. Also 
shown are the pseudo self-energies as well as self-energies of 
the baths. 



this end we employ a conserving approximation speci- 
fied by a generating Luttinger-Ward type functional <i> 
from which all self-energies are obtained as functional 
derivatives, = The building blocks of $ are 

the dressed Green's functions of pseudofermions G / (de- 
picted as dashed line), slave bosons Gf, (wiggly line), bath 
fermions Gc (solid line) and bath bosons Gh (curly line) 
and the vertices corresponding to hopping V, and ex- 
change interaction, /. 

In the strongly incoherent regime we are interested in, 
vertex corrections are not expected to change the be- 
havior in a qualitative way. They may, however, lead 
to quantitative changes. In this paper we would like 
to explore the leading behavior first, so that we may 
neglect vertex corrections for the moment. The lowest 
order terms of <I> in self-consistent perturbation theory 
in the bare coupling constants V and / are shown in 
Fig. The first one is the known generating functional 
of NCA, whereas the second one is new and involves the 
bath bosons. The corresponding self-energies are shown 
in Fig.^, for the pseudofermions (S/) and slave bosons 
(Ef,), as well as the bath fermions (Sc) and the bath 
bosons i^h)- We note that the impurity single-particle 
Green's function after projection (A oo) is related to 



(16) 



and likewise the impurity spin susceptibility is propor- 
tional to the bath boson self-energy 



(17) 



Cutting a pseudofermion line in each of the two diagrams 
of the generating functional. Fig. one finds two dia- 
grams for the pseudofermion self-energy 



as depicted in Fig. Likewise, the slave boson self- 
energy Efe, the fermion bath self-energy and the boson 
bath self-energy E/i are obtained by cutting the respec- 
tive Green's function lines in the two diagrams of $. The 
corresponding analytical expressions are given by 

^f^\ioj) = -V^TY,Gca{ii^')Gb{iLJ~iuj') , (19a) 

cr' ,Q 

xTj2Gha{t^)Gf„>{iuj + in) , 
n 

(19b) 

J^hii^) ^V^T^Gca{ii^')Gfa{in + iio') , (19c) 

a.Lu' 

^aaiiuj) ^ ~V^Tj2GfaiiuJ + i^)Gb{in) , (19d) 
o 

I];.„(zf7) = i/2^T,^,C, 

a, a' 

xTj2Gfa{iuj')Gf„,{tuj' + , 

(19e) 

where iu},iuj' and ifi are fermionic and bosonic Mat- 
subara frequencies, respectively; a, a' a = 1,2,3. 

Next one may transform the Matsubara frequency sums 
into frequency integrals along the branch-cuts of the 
Green's functions and perform the analytical continua- 
tion to the real frequency axis. The projection to the 
singly-occupied sector of Hilbert space may now be car- 
ried out. To this end the frequency arguments of the 
pseudoparticle Green's functions are shifted by the chem- 
ical potential A and the limit A — > cxd is taken. This yields 



(20a) 



(26), 



(18) 



(20b) 

Ef,(^ + iO)=^ / d^fiOAcAOGfA^ + ^ + iO)- 

(20c) 

where /(^) and n(^) are the Fermi and Bose functions, 
respectively, and Ac{^) and Dh{S,) are spectral functions 
of the fermionic and bosonic baths as defined in l|14|l . 

Since we incorporated the factors of and P into the 
definition, Ac and Dh are not normalized anymore, their 
total weight being given by V"^ and respectively. 
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The projected pseudoparticle Green's functions are ex- 
pressed in terms of their self-energies as 



Gf{uj + iO) = 



1 



a; + /i — Aq — S/(ti'-|-iO) ' 
1 

a; — Ao — Sb((jj + iO) ' 



(21a) 



where the (finite) energy shift Ao is determined by fixing 
the local charge 



= JdLue-f'^\j:^AfjLo)+A,{Lu)] =1. 



Here the subscript G specifies an expectation value in the 
grand canonical ensemble and Af(uj) = —^l'mGf{uj + iO), 
etc. 



The remaining self-energies Sc and E/j contain one 
pseudoparticle loop each and are therefore cx e~^^. The 
(21b) projected expectation value of any operator that vanishes 
in the Q = subspace is then given bjiSi 



{A) = hm ^ = lim e^'(A)c 



(23) 



(22) using 1(221). It follows that 



Ee,.(cj + iO) = y2 J ^^^-K yGf,{i + UJ + zO)Ab(0 - Af,{i)Gb{i - uj - iQ) . 
With the help of we find the imaginary part of the impurity Green's function in the compact form 



Im Gjmp,cr(t^ + iO) 



From 1)191;^ one finds after analytical continuation and projection 



d^e-f^^Af,{^ + uj)A,{£,). 



I],„(c. + zO) = - ^ r^„,r^,, j die-P^ [Af^iOGfa- {£. + + iQ) + G^^ - - iO)Af,,{0 

a. a' 

The impurity susceptibility is obtained from H17I26|) as 

Im X^mpA^ + ^0) = ^-'^"v / d^e-^^Af„{^ - u;)Af,, (O- 

(T.u' 

I 



(24) 



(25) 



(26) 



(27) 



Equations (f7|)- (|13|l . together with the "impurity solver", 
Eqs. (Pl l9abc)- H22l25l27|l have been solved self- 
consistently. Starting with given initial values of the 
fermionic and bosonic bath and pseudoparticle spectral 
functions, Ac{£) and Dh{£,), the first approximation to 
the pair of impurity Green's functions Gioc and xioc as 
well as the pseudoparticle spectral functions is deter- 
mined. Using the identities 



Gi, 



Xlo 



y I 



k 



^ ..-1 _ T2. 



(28) 
(29) 



that follow from equations H2I7I8I10|) and p7l9lll|l . the 

new bath spectral funtions Ac = —^ImV'^Gc and Dh — 
— ilm/^G/i may be deduced. With these and the up- 
dated pseudoparticle Green's functions one determines 
new GioctXiocG f ,Gb with the help of the impurity 
solver. The iteration is continued, until convergence is 
found to the desired level. This process is found to con- 



verge well in the temperature regime T > 0.04t using 
a nearest-neighbor tight-binding model, where t is the 
hopping amplitude. At lower T a solution could not be 
found any more. In the following we will present the 
results of the numerical evaluation, before discussing in 
detail the reasons for the breakdown of the solution in 
the low temperature domain. 



IV. RESULTS 

A. Local spectral function: pseudogap and 
non-Fermi liquid physics 

The most striking result of our work is the appear- 
ance of a pseudogap in the local electron spectral func- 
tion Aioc{'^) at small hole doping and low temperatures. 
Fig. 121 shows how the pseudogap starts to form when the 
exchange interaction J is switched on, for 5 = 0.015 and 
T = O.lt. In the limiting case of J = 0, corresponding 
to the Hubbard model in the limit U oo, Aioc{lo) is 
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FIG. 2: The local spectral function plotted versus frequency 
for four different J/t = 0, 0.1, 0.2 and 0.3 and T = O.lt for 
doping level of 5 = 0.015. The evolution of the pseudogap 
of width J is clearly visible. The zero of energy is set at the 
chemical potential ^. The inset shows temperature depen- 
dence of the local spectral function at the doping level 4% 
and for J = 0.3t. 



FIG. 4: The effective chemical potential fi^s vs doping for 
J = 0.3t and T = 0.06t (left scale). The dotted hne shows 
the non-interacting chemical potential fio- Open circles mark 
the estimation for the Kondo temperature Tk vs doping as 
calculated from Eq. 1311 (right scale). The arrow marks the 
position where Tk is equal to J. Only in the regime where 
Tk is larger than J, the solution shows the onset of a Fermi 
liquid phase. 
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FIG. 3: The local spectral function plotted versus frequency 
for T = 0.06t and J/t — 0.3 for various hole-doping concen- 
trations 5. The inset shows the characteristic temperature T* 
where the pseudogap opens (for the definition see the main 
text) . 



characterized by a broad maximum below the Fermi level 
{uj — 0) interpreted as the lower Hubbard band, and a 
narrow peak ( "quasiparticle peak") above uj = 0. As J 
is switched on, the quasiparticle peak disappears rapidly 
and the weight under it appears to be shifted a distance 
~ J below the Fermi level, forming a peak-dip-hump 
structure. The width of the pseudogap appears to scale 
with J. At the same time the spectral function develops 
a tail above to = reaching far (~ t) above the bare band 
edge. It is instructive to observe how the pseudogap dis- 
appears for a given J = 0.3i at T = 0.06i with increas- 
ing doping level (Fig. |3J). The pseudogap vanishes and 
the quasiparticle peak begins to appear at dopings above 
6 « 0.1. We note in passing that the bulk of the spec- 



tral weight in the lower Hubbard band is shifted rigidly 
with the chemical potential and only a section of width 
~ 4max( J, St) at the chemical potential is changing with 
the doping. 

The formation of the pseudogap at a low doping 6 = 
0.04 and fixed J = 0.3t as the temperature is lowered 
from T — 2J down to T = 0.2 J is shown in the inset 
of Fig. [3 In order to quantify the appearance of the 
pseudogap for given (5 as a function of T one may define 
the temperature T* at which the curvature of Aioc{i-^) at 
Lu = changes sign from negative to positive values as T 
is lowered. In the inset of Fig.Othe T* values determined 
in this way are plotted versus 6. T* is seen to drop rapidly 
with S, tending to zero a,t S 0.15. These results are 
reminiscent of what is seen in ARPES experiments on 
high Tc superconductorsii. 

How is the pseudogap generated? The clue to this 
question lies in the behavior of the effective chemical po- 
tential /ioff — fi — Rel](0), as a function of doping. In 
Fig. ^ ^off is shown at a low temperature T = 0.06t, 
in comparison with the bare chemical potential /io (of 
the noninteracting system). At doping levels 5 > 0.2 one 
finds that /icff coincides namely with fiQ , a necessary con- 
dition for Fermi liquid behavior. Upon lowering the dop- 
ing concentration, /XeS is seen to grow until at (5 « 0.02 
the upper edge of the bare band is reached (the zero of 
energy is fixed at the center of the tight- binding band). 
In fact /Zcff moves above the bare band, signalling the 
availability of states even above the latter. By contrast, 
the bare chemical potential remains located in the cen- 
ter of the band, approaching //q = in the limit S ^ 0. 
The fact that HeS is moving up towards the upper band 
edge for (5 ^ is a strong and unequivocal signal of non- 
Fermi liquid behavior - it is only possible for a highly 
incoherent metal with a large ImE. It is interesting to 



recall that in DMFT for the Hubbard model (which in 
the limit C/ ^ oo is identical to the t-J model for J ^ 0) 
one finds Fermi liquid behavior at low temperatures, and 
/.teff = Ho- Even at not so low temperatures (T > 0.06i) 
for J = /^cff follows /io except at rather low doping 
values 5 < 0.05, where a strong temperature dependence 
appears. 

Similar behavior has been found in Refs. for the 

Hubbard model at intermediate coupling. In Ref. a 
the dynamical cluster approximation involving up to 64 
sites was employed and the mean-field equations were 
solved by QMC simulation and the maximum entropy 
method, to effect the analytical continuation from imag- 
inary to real frequencies. Maier et al.^^ interpreted the 
pseudogap found in their spectra as generated by finite 
range antiferromagnetic correlations on the cluster or as 
RVB physics. Since the results are so similar to ours, and 
within our approach finite range AF correlations or the 
formation of intersite singlets are not included, we sug- 
gest that their pseudogap is created by the same mech- 
anism we identify as being responsible for our pseudo- 
gap: incoherent fluctuations (see above). Stanescu and 
Phillipsii used a two-site cluster approach to derive non- 
local DMFT equations. The quantum impurity model 
was solved by an adaptation of the non-crossing approx- 
imation. Again the results for the spectral functions are 
similar to ours. The authors claim that an effective low- 
energy model cannot be defined, as low and high-energy 
sectors are mixed in a dynamical way. We do not see 
any reason for such an unusual situation, neither from 
their paper nor from outside arguments. Rather, in the 
limit U 3> t, or more precisely, if U is strong enough 
to generate a Mott gap, the separation of the lower and 
upper Hubbard band is well defined, and a projection 
onto the lower band is justified. In Ref. [13 the ap- 
pearance of the pseudogap is attributed to short-range 
(nearest- neighbor) correlations, limiting the phase space 
for low-energy excitations. These correlations are iden- 
tified as orbital ring currents flowing between three ad- 
jacent sites. Since such effects are not included in our 
calculation, and we nonetheless find a pseudogap and a 
violation of Luttinger 's theorem, very similar to Ref. [l7j | , 
we conclude that the (and sketched above) interpretation 
given in Ref. fl^, is not correct. 

We conclude that the behavior found in our scheme 
for low doping, namely pseudogap and non-Fermi liquid 
physics, is a generic feature of an incoherent metal. We 
have found this incoherent state to be quite robust, e.g. 
against changes in band structure. It is worth mention- 
ing that ParcoUet and Georges^i recently studied & t-J 
model with random J, which is equivalent to our EDMFT 
equations for the Bethe lattice. They did not find indica- 
tions for a pseudogap. We believe the reason is that they 
employ slave boson mean-field theory, and thereby miss 
the incoherent part of the spectral function. A similar 
spin model has been considered before by Sachdev and 

At larger dopings the solution shows the onset of a 
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FIG. 5: Temperature evolution of local spectral function 
for doping level 5 = 0.24. In the main panel, the arrows 
show the position of the chemical potential while the inset 
shows spectra with the chemical potential fixed at tj = 0. 
The evoultion of broad quasiparticle peak above the Fermi 
level is clearly visible. 



Fermi liquid phase, which we now proceed to discuss. 
First we show in Fig. |5l the local spectral function at 
5 — 0.24 and J = 0.3i. With increasing temperature the 
quasiparticle peak broadens and the chemical potential 
shifts to lower energies. The value of Aiui) at the Fermi 
level increases with falling temperature and tends to a 
limiting value as T — + 0. 

At large doping the exchange interaction J is unimpor- 
tant and the EDMFT model reduces to an Anderson im- 
purity model. We may estimate the hybridization width 
F of this model from the density of states of the fcrmionic 
bath at the Fermi level (w = 0): 



F = ttA,{uj = 0). 



(30) 



The energy of the local orbitals Ed, according to ©. is 
given by the chemical potential Ed = — /i. An estimate of 
the Kondo temperature is obtained from the conventional 
expression Tk = V DT exp{^^^^) as 



Tk = ^DnA,{0)exp{- 



2Ac(0)' 



(31) 



Fig.0]shows Tx/t as a function of S for the low tempera- 
ture T = OAt, using D — 2t. The Kondo temperature is 
seen to fall strongly with decreasing 5 even at the highest 
value i5 = 0.24, and approaches zero rapidly in the pseu- 
dogap regime. The value where Tk = J is indicated. In 
the regime Tk ^ J one expects the exchange interaction 
to be of increasing importance, such that the interpre- 
tation in terms of an Anderson impurity model loses its 
meaning. 

In the Fermi liquid regime the imaginary part of the 
self-energy of G/oc is expected to vary as 



Im T,{u - iO) - < LU^ + {nTf /T\ 



n2 

'K 



(32) 
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FIG. 6: Quasiparticle renormalization amplitude Z plotted 
vs temperature for various doping concentrations. The inset 
shows the imaginary part of the self-energy at zero frequency 
as a function of temperature. 



where the Kondo temperature Tk plays the role of the 
renormalized Fermi energy. The quadratic dependence 
is expected to hold for u;,r ^ Tr:. The inset of Fig. 
shows Im S(0) as a function of T for doping levels from 
(5 = 0.3 down to 0.01. A limiting quadratic tempera- 
ture dependence is not seen since the lowest temperature 
reached in our evaluation is above Tk (or, for b > 0.2, 
only slightly below Tk)- 

However, for S — 0.24 and 0.3 behavior consistent with 
Fermi liquid theory would smoothly match the results 
shown; For smaller doping, in particular around S ~ 0.1, 
Im S at T ~ 0.03^ is so large that it is impossible to 
connect this behavior smoothly with a Fermi liquid be- 
havior below Tk — O.lt. At still smaller 6 Im E is 
seen to decrease with doping, due to the formation of the 
pseudogap. 

As a further indication of Fermi liquid behavior we 
evaluate the quasiparticle weight factor Z defined as 



1 



aReEx-i 



(33) 



Fig. shows Z as a function of temperature for S ~ 
0.02—0.3. A finite quasiparticle weight in the limit T — > 
would signal Fermi liquid behavior. It is seen that only 
for the highest doping levels S = 0.3 and 0.24 would an 
extrapolation to T = give a finite value. For smaller 
values of 6 the Z-factor appears to drop rapidly with 
decreasing temperature, possibly extrapolating to zero. 



B. Pseudoparticle spectral functions 

The pseudofermion spectral function Af{uj) a,t S = 
0.015 and T = O.lt is plotted versus uj/t in Fig. □ for 
values of J/t from to 0.3. While at J = A/ (w) is 
characterized by a narrow peak at a frequency lu T and 
of width ~ T, increasing J leads to a rapid asymmetric 



FIG. 7: The pseudo-fermion spectral function plotted vs fre- 
quency for four different values of J. The inset shows the evo- 
lution of spectra by doping the system at constant J — 0.3t. 
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FIG. 8: The pseudo-boson spectral function for the same 
parameters as used in Fig. [7| 



broadening of the peak, of width Aw ~ J. Although in 
the limit T ^ for general reasons one expects Aj(uj) 
(and also Ah[iS)) to acquire power-law divergent behavior 
at the threshold w = 02&, the temperature T = O.li is too 
high to show the asymptotic behavior. At large doping, 
6 > 0.2, when the Kondo temperature Tk as defined in 
(|31|l is larger than J, Af{uj) is hardly affected by J. The 
doping dependence of A/(w) at J = 0.3t, as shown in the 
inset of Fig. [7| is weak. The characteristic energy scale 
is max( J, Tk) ~ J up to the highest doping of (5 = 0.24, 
and hence is independent of 6. 

The pseudoboson spectral function shown in Fig. |S1 is 
roughly speaking a mirror image of the lower Hubbard 
band. As J is switched on spectral weight is pushed from 
below the threshold at = and from the far end of the 
Hubbard band into a peak at a; ~ J, emulating the peak- 
dip- hump structure in Aiodi^) in the pseudogap regime. 

Both in the case of J = and for S > 0.2 a sharp 
quasiparticle peak is observed to form in Ah^uj) at w = 0. 
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FIG. 9: The local dynamical spin susceptibility plotted vs 
frequency for four different J/t and doping concentration S — 
0.01. 



FIG. 10: The spin self-energy M(u}) plotted vs frequency for 
four different J/t. The horizontal dotted lines mark the value 
zj, where z = 4 is the coordination number. 



The peak is suppressed at temperatures T 3> Tk, which 
is why it is not seen in Fig. |S1 In contrast to Af{uj), 
Ai,{lu) is strongly doping dependent, as shown in Fig.|Sl 
At the moderately low temperature oi T — O.lt and for 
J = 0.3t, the pseudogap feature at small 6 is seen to 
smoothly cross over to the quasiparticle peak at large 
doping. 



C. Dynamical spin susceptibility 

The dynamical spin susceptibility is expected to re- 
veal how the character of spin fluctuations dependends 
on doping and exchange coupling constant J. In Fig. |3 
the imaginary part of x/oc(w)/w is shown at low doping, 
6 = 0.01, and low temperature, T = 0.09^, for various 
values of J/t ranging from to 0.3. As J is increased, 
the peak of I'm{uj)x/LL! broadens and the width is seen to 
be given by Aw ~ J. The real part Rex{oj = 0) = x'(0) 
decreases with J as shown in Fig. O However, there 
is no trace of a pseudogap in Imxioci^)- The pseu- 
dogap reveals itself in the spectrum of the self-energy 
M{oj) of magnetic excitations, as shown in Fig. 1101 where 
lmM{uj) = M"{uj) is observed to develop a gap for 
Lo < J . As analyzed in Section ^ the pseudogap is 
caused by large values of x'(^)i which force a redistri- 
bution of spectral weight in M"{Ld) by way of the self- 
consistent feedback of x'(<^) into M"{u)). In Fig. ^] the 
momentum- resolved spin excitation spectrum, x"^(lo)/lo 
is shown for J = 0.3f, 5 = 0.01 and T O.lt. Whereas 
a pronounced gap exists at g-values away from the anti- 
ferromagnetic wavevector Q — (tt, tt), near Q the gap is 
filled in. This is due to the fact that in the region of q- 
space around Q not only M" {lo) is small for w < J, but 
also the real part of the denominator of Xqi'^) vanishes, as 
M'(w)-|- Jg — > for q-^ Q and as the transition to the an- 
tiferromagnetically ordered state is approached. Conse- 
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FIG. 11: The momentum dependent spin susceptibility along 
the (0, 0) — (tt, tt) axis plotted as a function of frequency for 
doping concentration 5 = 0.01 and temperature T — O.lt. 



quently, the ratio M"{uj)/\xq{uj)\'^ develops pronounced 
peaks at |a;| J rather than a pseudogap. In the lo- 
cal susceptibility the contribution from q k Q tends to 
fill in the pseudogap, which is therefore not discernible in 
Fig. El The effect of approaching the ordered state is also 
observed in the real part of Xg(^)j shown in Fig.^] The 
static q-dependent susceptibility Xg(0) is seen to grow by 
two orders of magnitude as q is varied from g = to 
q — Q. This behavior reflects the effect of a large spin 
correlation length ^, defined through 



zjr' + (q-Q)'' 



(34) 



for q « Q (z = 4 is the coordination number and length 
is measured in units of the lattice constant). In Fig. 1121 
the inverse correlation length is plotted versus T/t for 
J = 0.3t and for various doping levels. For comparison, 
the theoretical result for the Heisenberg model (two loop 
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FIG. 12: The inverse of dynamic correlation length plotted 
as a function of temperature for various doping levels. The 
curve for 5 = is taken from Ref. [2^ and corresponds to the 
two-dimensional Heisenberg model. 



posed self-consistency of the EDMFT equations has led 
to a drastic renormalization of the structureless tight- 
binding density of states. In fact Ac(a;) reflects the 
structure seen in Aiocii^') to a large degree: On the one 
hand, the quasiparticle peak at large doping and on the 
other the pseudogap at small b. For comparison we show 
Aioc(w) in Fig. El as well. 

A similar trend is seen in the case of the spectral func- 
tion of the bosonic bath Dh{uj), as is apparent from 
Fig. 1141 Here we also see from the comparison with 
Xioci^) ^ large degree of similarity. 

The total weight under the spectral functions Ac{uj) 
and Dh^Lo) is equal to the squares of the coupling con- 
stants V"^ and respectively. As shown in Appendix 
IbI the coupling constant is fixed by sum rules and is 
given by 



-S). 



(35) 
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FIG. 13: The fermionic bath spectral function Ac for two 
different doping levels S = 0.01 and 5 = 0.18 at J = 0.3t 
and T = O.lt. The local spectral function is also shown for 
comparison. 



In contrast, a similarly simple relation does not hold for 
P. However, may be related to xioc and M as follows: 



duj 
7^ 



Im[M(L^-iO)-xroc(^-*0)] ■ (36) 



It turns out that the numerical evaluation yields 



/2«2j2(l-,5). 



(37) 



The first moment of the eigenfrequencies ojq of the 
bosonic bath, is given by the f-sum rule. 



-E-. = |^/^V"^"o». (38) 



order RG of the nonlinear sigma model) given in Ref. 
hM (limit 6 = 0) is shown as well. It appears to connect 
smoothly to the curve for 6 = 0.02. Fig. E|also serves 
to show that the numerical solution ceases to exist at 
Xq(0) > 10^, as will be discussed in Section Ivl 



where ^e^) = / dee^Nj{e) and Nj{e) is the density of 
states (DOS) of Jq. 



D. Spectral functions of the fermionic and bosonic 
baths 

The spectral function Ac{uj) of the fermionic bath is 
shown for J = OM and at T = O.lt in Fig. US) The im- 



E. Thermodynamic properties 



The thermodynamic potential fl within EDMFT can 
be expressed in terms of the impurity free energy Qimp 
and contributions from the fermionic and bosonic baths: 



12 



iiu I k.cr q,Q J 



Performing the analytical continuation from imaginary frequencies to the real axis and expressing the momentum 
summations as energy integrals, (|39|l may be written as 




co/t 



FIG. 14: The bosonic bath spectral function for two dif- 
ferent doping levels S = 0.01 and 5 — 0.18 at J = 0.3t and 
T = O.lt. For comaparison, the local dynamic spin suscepti- 
bility is also shown. 



FIG. 15: Entropy per site as a function of doping 5 at various 
temperatures. Exact diagonalization results^** for the same 
temperatures are denoted by dotted lines while the open cir- 
cles correspond to the experimental data-^ on La2-a;Sri:Cu04. 



The impurity free energy is given by the shift of the 
chemical potential^, Aq, defined by 



Ao. 



(41) 



The entropy S — —(^^^ as a function of doping con- 
centration (5 for various temperatures is shown in Fig. 1151 
Even at the low temperature, T — O.lt, S is seen to be 
rather large ('^ 0.5 In 2), an indication for strong corre- 
lations and a rather incoherent state. The entropy of 
a noninteracting system at the same density would be 
about an order of magnitude smaller. The overall mag- 
nitude of S compares well with both the results of ex- 
act diagonalization^ for a small system and experimen- 
tal data for La2_2;Sr2;Cu04 29]. The calculated entropy 
shares the trend that it is reduced both at large doping, 
when the system crosses over to a Fermi liquid, and at 
smaller doping in the pseudogap phase. The quenching of 
the magnetic fluctuations by the incipient magnetic order 
as the antiferromagnetic Mott insulator is approached for 
(5 is qualitatively reproduced (note that for J = 0, 
S increases as 5 — > 0, and this behavior is obtained in 
DMFT calculations of the Hubbard model). 
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FIG. 16: Specific heat coefficient vs temperature for various 
doping concentrations. In the right panel we show results 
obtained by the exact diagonalization^. 



In Fig.^]the specific heat divided by the temperature 
is plotted versus T for various dopings 6 (left panel). 
For comparison the results of ED arc shown on the right 
panel. 
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FIG. 17: Solid line: derivative of the thermodynamic po- 
tential with respect to chemical potential 1 -I- (9f2/9/i)T (or 
equivalently doping vs. chemical potential). The contribu- 
tions from three different parts of the thermodynamic po- 
tential: impurity, electron Green's function (second term in 
Eq. 1401 and spin susceptibility part (last term in Eq. I4UII are 
shown separately. 




FIG. 18: The weigting functions for the two di- 
mensional square lattice can be expressed by ele- 
mentary functions as ^xx{x) = ^ [2|a;|i5(l — 
+2K{l-l/x^) -2Il{l-l/\x\,l-l/x^)] and <S>xy{x) = 

2 (f )^ [x^E{l - Xjx^) -K(\ - l/a;2))]Sign(a;). Here, K(x), 
E{x) and n(a;) are complete elliptic integrals of the first, 
second and third kind and x = e/(4f). 



Another thermodynamic quantity of interest is the par- 
ticle density n, given by 



(42) 



In Fig. El the doping & is plotted versus /i at T = O.li. 
As expected, 5 varies monotonically with /i, with positive 
curvature. 

The particle density may also be obtained from the lo- 
cal Green's function as n = 2Gioc,a{T = 0+). The result- 
ing values of n are indistinguishable from those calculated 
by differentiating fl, which provides a check for numerical 
accuracy within our conserving approximation. 



F. Transport properties 

The calculation of transport properties in EDMFT is 
facilitated by the observation that a momentum indepen- 
dent self-energy leads to a local current vertex function 
(in other words, the non-local parts vanish in the limit 
dimension d oo)^2*2ii2^. The optical conductivity is 
therefore given by the single-particle Green's function as 



iu' kcr 



(43) 

where = 2t sin is the bare current vertex. Using 
the fact that Gk depends on k only through Ck [see 
Q], and performing the analytical continuation to the 



real frequency axis one finds 



X duj 



A{e,uj')A{€,uj' + u;) (,44) 



where 



(45) 



and A{ek,ijj) = ilmGk(w — iS). 

Similarly, the off-diagonal or Hall conductivity in the 
presence of a magnetic field B perpendicular to the plane 
takes the formS*Mi2^ 



47r2e3 



^ xy — 

where 

and 
det(k) 



A{e,u;) 



=I]det(k)5(e-ek) 



dka 



af3 



(46) 
(47) 



dkadkp 
(48) 



The weight factors ^xx and are shown in Fig. E| 
One observes that for the simple 2D tight-binding lat- 
tice, ^xx is even function of energy while ^^y is an odd 
function of energy. 

The Hall coefficient TZh is defined as 



<^lxB- 



(49) 
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FIG. 19: T-dependence of the resistivity multiplied by doping 
S. The linear T behavior for high T flattens for 5 > 0.1 at 
a temperature of the order of J. For S < 0.1 the resistivity 
drops in the regime where a pseudogap opens. 



FIG. 20; T-dependence of Rh for J ^ 0.3t. For small doping 
and T — > 0, Rh approaches the value l/(|e|(5) expected for a 
single hole in a t-J model. 



For orientation it is useful to discuss the limit of low 
temperatures, assuming ImS](w) at uj ^ 0, and 
A(e, Lo) sharply peaked as a function of e at e — uj + fics- 
One may then do the integrations on e and uj in H44I) and 
(gni), yielding 



and 



n 



H 



|ImS(0)| 



1 ^xyi^-cs) 

2e ['^xx{^Ji<,s)Y 



(50) 



(51) 



with e = — |e|. We observe that TZh does not depend on 
ImE in this limit. In the Fermi liquid regime /^off — fJ-o < 
0, and consequently ^xyifJ-o) > Oi leading to a negative 
Rh < 0. 

By contrast, in the incoherent regime of the t-J model 
^cS is found to be positive, approaching the upper band 
edge for (5 — *■ (see Fig.EJ. Since ^xy{e) is negative for 
positive e, and TZh is seen to be positive (hole-like). For 
6 > 0.17, fj,cs changes sign and TZh turns negative. For 
the nearest neighbor tight-binding model, and assuming 
a linear variation of fXcS with S, fi^s = 4t(l — CS), TZh 
takes the simple form 



TZh 



n 1 
20^6 



; C > 0. 



(52) 



Using in addition the result for a single hole in the half- 
filled band2&, TZh ~ Vkl^' o^^^ finds by comparison C = 
^. For the conductivity one obtains in a similar way 



t6 



|lml](0)| 



0. 



(53) 



Although (|52|l and H53|) are in qualitative agreement with 
our numerical results, we emphasize that the assumption 



of small ImE(O) is not justified in the incoherent regime. 
A large ImE is actually necessary to obtain a /ies close 
to the band edge and therefore a positive sign of Rh- 

We now present the numerical results. In Fig. El the 
scaled resistivity PxxS/po, where po = fi/e^, is plotted 
versus temperature for values of S ranging from 0.01 to 
0.23. The curves form a narrow band meaning that the 
scaling p^x oc 1/S shown in (|53|l holds approximately 
(and ImE(O) is a weak function of 6). The values of 
the resistivity are rather high. In the pseudogap regime 
{S <C 0.1) the resistivity tends to turn downward for de- 
creasing temperature. By contrast, at higher dopings an 
upward curvature is observed, leading to a plateau at low 
T, before p begins to drop to lower values at still lower T. 
We cannot exclude that the plateau is an artifact of the 
NCA approximation. The Hall coefhcient is plotted in 
Fig. 1201 versus temperature, for values of 6 ranging from 
0.01 to 0.23. For small doping S < 0.16, TZh is always 
positive, approaching the expected valued l/(|e|(5) in the 
limit (5 ^ 0,T ^ 0. For doping levels S > 0.16 TZh is 
negative at low T, consistent with (|51() . and changes sign 
at higher T, similar to what is observed in experiment^. 



V. INSTABILITY OF THE EDMFT SOLUTION 
DUE TO CRITICAL FLUCTUATIONS IN d = 2 



We will now investigate the question of why no solu- 
tions of the EDMFT equations exist for low temperatures 
and small dopings. Within the EDMFT of the t-J model 
long-range antiferromagnetic fluctuations are not taken 
into account in a proper way. As a consequence the local 
spin excitation spectrum Xi'oci^) ~ IniX/oc(^^ — ^0) keeps 
a simple Lorentzian-type shape. On the other hand the 
static local susceptibility x!oc(0) i^^ dimensions) di- 
verges as In ^ when the transition to the antiferromagneti- 
cally ordered state is approached and the spin correlation 
length ^ — > oo. This in turn forces the slope of X/oc(^) 
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in the limit — > to diverge as In^ as well. Within 
the effective impurity model of EDMFT a steep slope of 
Xioci^) entails a large maximum of x'Lci^) ^max ^ J^ 
of value x'loci'^rnax) ~ x'loci^) ■ ^^^^ shown below, a 
maximum value of xHoA^rnax) larger than some critical 
value X/oc crit ~ 7 ' where the constant c depends on the 
density of states Nj{e) (see (|54|l '). leads to an unphysical 
pole in Xq(cj) at w = ujmax and q = Qmax- This in turn 
forces ImM(a; — iO) to change sign into an unphysical 
branch of the complex frequency plane. This is the point 
when a stable numerical solution cannot be found any 
longer. 

To demonstrate this behavior explicitly we consider 
now a flat density of states of spin excitations, 

Njie) = 5] <5(e - Jq) = i^0(4 J - \e\). (54) 



where the bandwidth 8 J has been chosen to agree with 
that of the tight-binding model, Jq = 2J(cosgx + cos(7j^). 
The local susceptibility as defined by may then be 
expressed analytically in terms of the self-energy M{lu) 

1 , AJ + M(uj-iO) , „ , , 
^'-(--Q)-8J^'- -4J-^M(.-.0) -^ +-^ - 

Inverting this relation one finds 

v + l |i;P-l-2ro" , , 
M{uj - iQ) = 4 J = U— , (56) 



w-1 



where v = exp(8 Jx) 
V, given by 



The imaginary part of 



v" = exp(8 Jx') sin(8 Jx") 



(57) 



will change sign as x"{^) increases with increasing w, 
if 8Jx" > TT. By (|56|) . this will lead to a sign change 
of M"{u) — iO) from negative (stable) to positive values. 
How can x" and M" both be positive? This is possible 
since Xq{^ ~ *0) develops a pole in the physical domain, 
— I Jql < 4 J, at finite w — ujmax, giving a contribution to 
Xioc with the "right" sign. The instability occurs at finite 
frequency and thus is not easily interpreted as a physical 
phenomenon. 

In the numerical treatment we found that a convergent 
solution cannot be obtained when the stability criterion 



// 

Xlo( 



< 



J 



(58) 



is violated. The constant c takes the value tt/8 for the 
fiat DOS and a value ~ 0.3 for the tight-binding model. 

We emphasize that this instability is not an artifact 
of the method of solution of the impurity model but is 
a generic feature of the EDMFT equations in two di- 
mensions. We argue that whenever the ground-state at 
T = is ordered, the self-consistency scheme has to break 
down below some finite temperature. This argument is 



not only relevant for our calculation but should be rel- 
evant for other applications of EDMFT which have fo- 
cused on discussing the possibility of novel quantum crit- 
ical points in the presence of two dimensional magnetic 
fiuctuations'^^-'^'*. While our reasoning does not apply di- 
rectly to the quantum critical point, it strongly suggests 
that no solution exists on the ordered side of the phase 
diagram, casting doubt on the applicability of EDMFT 
also at the quantum phase transition. Our formal argu- 
ment starts from the observation that in two dimensions 
no phase transition (of first or second-order) is possible 
for T > within EDMFT, since in a hypothetical or- 
dered phase the local susceptibility would diverge due 
to the presence of Goldstone modes - in this respect, 
the EDMFT approach obeys the Mermin- Wagner theo- 
rem. We mention in passing that even in the case of 
Ising symmetry a second-order phase transition is not 
possible, as within EDMFT the longitudinal fluctuations 
would diverge at the critical point; however, in this case a 
first-order transition towards an ordered phase for T > 
cannot be excluded on general grounds. Indeed, a first 
order transition has been found by Sun and Kotliar— 
and Zhu, Grempel and Si^^ for an Ising-coupled Ander- 
son lattice or Kondo lattice, respectively. Assuming that 
for T = the system is magnetically ordered, the lo- 
cal susceptibility Xioc(O) ~ In^ will grow steadily as T 
is lowered where ^ is exponentially large, ^ ~ e^^ , for 
T <^ E* and E* can crudely be identified with the mean- 
field transition temperature. However, we have shown 
that within EDMFT x"(^) is bounded from above by 
the requirements of self-consistency. How can this be 
reconciled with large x;oc(0) = / l212LlMd^ _ in^? The 
only possibility consistent with the Kramers-Kronig re- 
lation is that x"(^) is constant down to an exponentially 
small energy scale E'^'^ ~ where z is some positive 

exponent. For sufficiently small T, E'^'^ will be exponen- 
tially smaller than T. At this point we have to ask the 
question whether the solution of the effective impurity 
model can result in a scale exponentially smaller than 
T. We think that this is extremely unlikely and conclude 
therefore that no solution can exist for sufficiently small 
T, consistent with our results and also with QMC simu- 
lations by Burdin et al.'^^ of a model equivalent to ours 
in the limit of zero doping. 

Eq. also shows how the pseudogap in AI{uj) 

emerges from the self-consistency of x and M. The ab- 
sorptive part of the self-energy M, as seen from (|57|) . is 
exponentially small in the regime where 



8Jx'(c^) » 1. 



(59) 



From the numerical results in Fig. 1101 one sees that H59() 
is satisfied if < cmJ, where cm is a constant of 
order unity, which depends on the DOS Nj{e). Thus, 
the pseudogap is found to develop as a consequence of 
the increase of x'(0) jln^ with growing ^, in two- 
dimensions. We stress that a relation similar to (|56|l be- 
tween M{u!) and x('^) holds whenever the DOS Nj{e) 
is finite at the band edges, which is a signature of two 
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dimensions. In this case the conclusions drawn above 
remain valid when -!j is replaced by ANj, the DOS 
jump and x is replaced by x — Xreg, where Xregi^) = 
J de{Nj{e) - A7Vj)/[Jq + AI{lo)]. 

VI. CONCLUSION 

The physics of the doped Mott-Hubbard insulator is 
governed by the interplay of the motion of holes and the 
antiferromagnetic fluctuations of the spin background. In 
this paper we have used a local approximation scheme to 
describe both the constrained hopping of holes and the 
quantum spin fluctuations in the paramagnetic phase on 
an equal footing. The local approximation becomes ex- 
act in the limit of infinite coordination number of the 
underlying lattice and is known as Extended Dynami- 
cal Mean Field Theory. Rather than studying the model 
in this limit, we take the point of view that in finite di- 
mensions the approximation of neglecting the momentum 
dependence of the single-particle self-energy and the J- 
irreducible spin susceptibility may still be useful. Here we 
have applied this scheme to the two-dimensional square 
lattice with nearest-neighbor hopping and exchange in- 
teraction. We expect that the approximation should 
work in a regime of temperatures and doping concentra- 
tions where incoherent fiuctuations dominate and wash 
out any of the collective effects sensitive to the system 
dimension, such as long-range antiferromagnetic order or 
superconductivity. 

In the regime of temperatures above T ^ O.IJ and for 
doping levels of 0.01 < S < 0.3, we indeed find a highly 
incoherent phase, with a broad distribution of spin ex- 
citation energies, a high entropy and a large electrical 
resistance. Most strikingly, the local single-particle spec- 
tral function, which is characterized by a narrow peak 
above the chemical potential for i5 > 0.25, develops a 
pseudogap as 6 is reduced down to the few percent range. 
The appearance of the pseudogap is related to a dramatic 
shift of the effective chemical potential from its noninter- 
acting (i.e. Fermi liquid) value near the center of the 
lower Hubbard band to the upper band edge. The shift 
persists down to the lowest accessible temperatures of 
T « O.IJ and constitutes an unequivocal signal of non- 
Fermi liquid behavior in the regime 0.01 < (5 < 0.2. The 
single-particle pseudogap is accompanied by a gap in the 
spin excitation spectrum for momenta not too close to 
the ordering wave vector Q = (tt, tt). 

The Hall transport is found to be hole-like, the Hall 
constant tending to large positive values cx 1/S as the 
doping is reduced. At large dopings and low tempera- 
tures Fermi liquid type behavior is recovered. 

These results are encouraging and give rise to the ex- 
pectation that the present EDMFT scheme is able to 
capture the main features of the t-J model in the inco- 
herent regime. At lower temperatures and small dopings 
one should expect the closeness to the antiferromagnetic 



transition at T = and d < 6c to play an important 
role. We indeed find that the EDMFT equations stop 
having a physical solution below a limiting temperature 
of T « O.IJ. We are able to trace this behavior to an 
intrinsic lack of structure in the spin structure factor of 
the effective impurity model, which is ultimately due to 
the insufficient treatment of critical fluctuations in the 
EDMFT model. It is likely that similar limitations apply 
to other applications of the EDMFT in low dimensional 
systems. 

In conclusion, we emphasize that within the present 
local approximation scheme neither effects of finite- 
range, slowly fluctuating antiferromagnetic or supercon- 
ducting domains, nor local singlet formation or similar 
short-range correlations are included. Nonetheless, the 
strongly incoherent fluctuations characteristic of our ap- 
proach (in this case of the spins, but one could imagine 
similar effects e.g. in the superconducting sector) suffice 
to drive pseudogap formation, a violation of Luttinger's 
theorem and a hole-type Fermi surface in the proximity 
of a Mott insulator. 
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APPENDIX A: EDMFT DERIVATION 

In this appendix we derive the EDMFT self-consistent 
equations for the t-J model using the cavity method. 

To treat the no-double occupancy constraint of the t-J 
model, we will add a local Coulomb repulsion term ex- 
plicitly and take the limit U —f oo at the end. In this 
approach, the electron creation (destruction) operators ct 
(c|) obey the usual fermion anticommutation relations. 
The resulting Hamiltonian is the so-called extended Hub- 
bard model 

ij.a i ij 

(Al) 

It is straightforward to extend the theory to other non- 
local interactions like non-local Coulomb repulsion, but 
since we are mainly interested in the effect of magnetic 
fiuctuations we will neglect other terms in the Hamilto- 
nian. 

For simplicity, let us assume there is no long-range 
order (i.e. the system is in the paramagnetic state). Let 
us start the derivation of the EDMFT equations with the 
action corresponding to the Hamiltonian (|A1I) : 
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S 



dr 



(A2) 



The action can be divided into three parts: the on-site part for the chosen site {So) 



So ^ dr 



(A3) 



the inter-site interaction between the chosen site o and the rest of the system (AS') 

AS = dr y2^''^io4cri'^)(^ocr{T) - Ci4cr (t) + - ( Jio + Joi) Si(r) • So(t) 

Jo , ^ ^ 



(A4) 



and the lattice action in the presence of the cavity (5^°-'), which is equal to the original action HA2|I with site o 
excluded from all summations. 

The series expansion in the coupling between the central site and the rest of the system can be expressed as 

Z = y Dci^Bcoa J n DcIDc,^ exp (^-So - 5^") - AC{T)dT^ (A5) 
= J DclDcoaeM-So)Z^"^ (^l - j\ACiT))^°UT+ ^ dn dT2(r,A/:(ri)A£(r2))(°) + ... j , (A6) 

where AS" = A£(r)dT and ()(") means the average over the cavity action S^^^ . In the second line we have integrated 
out all fermions except for site o. 

The first term linear in AC vanishes, since the average of each spin (Si(T)) = is zero by the assumption of no long 
range order in the system. For the broken-symmetry phase, the spin operator has to be replaced with its deviation 
from the average value S.^ Si — (S.^) and the derivation can proceed along the same lines. The second term in the 
series expansion reads 



i r dT2{TrACiTi)AC{T2))^"^ = 

^- Jo Jo 

^tio4ai'^i)^oaiTi) -I" ^oj^ct ('^1 )Cicr (^1 ) ^ ^ Joi So(ri) • Si(Ti) 



/3 ./3 / 

dTi / dT2 ( Tr 

Jo \ 



'^tiocl^{T2)Coa{T2) + ioi4<T ("^2 )Cjo- (t2 ) - Jio Sj(r2) • So{t2) 



(0) 



(A7) 



It is crucial to observe that there is no interference between the kinetic and the spin term since the average of the 
correlation function (c,;(j(ri)Sj(T2))*'°'' vanishes. The leading-order term in the effective action thus reads 



SeS = So- / / dridr: 



(0) 



clA'Tl)^Uotoj iTrCia{Ti)c^^^{T2)^ Cqo- (r2 ) + So(ri ) ^ ^ JioJo, (T^-Si (ti )Sj (r2 )) So{t2) 

m 



Within EDMFT both terms are equally important and apart (in Manhattan distance) giving a contribution of 
are of order 1 in the 1/d expansion. The two-point order 1/d. The prefactor or is proportional to l/d. 
Green's function and the susceptibility scale as l/d'*"^'/^ while the double sum gives d^ and the net result is there- 
since t and J fall off as l/\/d. Furthermore i and j are fore of order 1. 
neighbors of site o and are thus at least 2 lattice sites 

Further it follows from the Linked Cluster Theorem 
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that only connected n-point correlation functions appear 
in higher-order terms of the effective action. Since they 
have the usual dependence on 1/d, all but the first term 
vanish in the limit d oo. For instance, the next-order 
term would involve 3-point connected correlation func- 
tion Xijk ^ (S^SjS^) or C'ijk ^ I^Sfc^jCk^ that scale like 
l/(fl»-j|/2g;|i-fe|/2^ When all three variables i, j and k are 
different, the correlation function is of order 1/d^ since 
all three sites are neighbors of o. The prefactor or Jt^ 
is proportional to l/d^/^ while the sums give d'^ . The 
term is thus of order 1 / ^/d. If i = j but distinct from k 
the correlation function is of order 1/d while sums give 
d^ and the net result is again of order 1/^/d. Higher- 
order terms fall off faster than l/\/rf. Thus, in the limit 
of large d all but the first term ljA8|l can be neglected and 
the effective action becomes 



ScB^ U no]iT)noi{T) 







/3 rP 

dTi dT2cl^{Ti)g^^{Ti - T2)Coa{T2) 

Jo 

^^n^ rfT2So(Ti)x-'(ri-T2)So(r2) (A9) 



where 



Xo 



(AlO) 



The Weiss fields are thus determined by the cavity 
Green's function G^°'' and the cavity susceptibility Xij'' ■ 



The absence of interference between the kinetic and spin 
terms in (|A8|I also leads to separate equations for both 
cavity quantities 



a 



(0) 



(0) _ _ -1 

Xij — Xij XioXoo Xoj- 



(All) 



With the power-counting arguments one can show^'^ that 
in the limit d — > oo and EDMFT scaling the single- 
particle self-energy Y^{iuj) as well as the double particle 
self-energy M{iuj) become local quantities, i.e.. 



X^(i^) 



1 



1 



,/q + M {iuj) ' 



(A12) 



Inserting the definitions (|A12(I into (|A11|I and combining 
with (|A10|) we finally obtain the self-consistent conditions 



Qo' = s + 



Xioc- (A13) 
These relate the Weiss fields to the local quantities com- 
putable from the local action (|A9(I . The system of equa- 
tions is thus closed. 

For practical computation, however, it is convenient to 
have a Hamiltonian representation of the local effective 
action HA9|I . Since it includes retardation effects through 
frequency dependent Weiss fields, it is necessary to intro- 
duce auxiliary degrees of freedom describing the baths. 
The one-particle character of the Weiss field Q^^ can be 
represented with the fermionic bath while the two parti- 
cle field has a bosonic nature and dictates bosonic 
bath. One of the possible choices is 



H = Y^ Ekcl^Cka + VY^{cl^Coa + C^Cfe^) - ^ ^i cl^Coa + Uno^UQ^ + ^^'^J • hq + / ^ • (h, + hLg), (A14) 



kij 



where corresponds to a vector-bosonic bath with the commutation relations [/i^, h^}] — Sqq'Sap- The corresponding 
action 

f/3 



^-^o + ^ drY. cUr)i^+Ek)ck.+Vcl{T)coAr) + Vcl{T)ckAr) 



drj: 



ht(r)(— + LOq)\{T) + /h,(r) . S<,(r) + /S<,(r) • hl^ir) 



(A15) 



is quadratic in Cka and h^, and therefore both baths can be eliminated leading to 



S^So- f dr, r dr2j2cl{r,)iJ2v'^^)coAT2)~ f dr, f dT2So{T,){Y,l'4^)^o{r2). (A16) 
^0 ^0 ^ ^. + Jo Jo „ d7+^i 



r 



This action is identical to the effective action (jA9p pro- vided that the following relations hold 

'A 



fe dr 



Ek 
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or equivalently 

Finally, combining (|A13p and l|A18(l yields 



G 



loc 



Xlo 



iLu + fi-j:-v^Gc , 

M + I^Gh , 



(A18) 



(A19) 



which coincide with H10|) and (|ll|l . 



APPENDIX B: SUM RULE CONSTRAINTS ON 

V'^ AND J2 

Within EDMFT, the couphng parameters V and / 
defined by lO, describing the hopping on to, and the 
exchange interaction with the impurity, are determined 
self-consistently. Interestingly, sum rules completely de- 
termine the values of V and partially constrain the values 
of /. We start with the single-particle hopping V. Defin- 
ing a complex variable z = uj + fi + S(a; — iO) one may 
write the EDMFT self-consistency condition (7), using 
(2), (8) and (10), as 



H{z) 



f D{e)de 



1 



z-V^Gcitu-iO) ' 



(Bl) 



with D{e) the DOS of the tight-binding band ek and 
Gc{(^) the fermionic bath Green's function defined in 
(12). Solving JbTJ for V'^Gc and taking the hmit Cij ^ oo 
one finds 

lim [luV'^Gc{uj - iO)] = V"^ = lim^^^oo i^iz - -nTzr) 

= lim„^eo7(e'). (B2) 

Here the zero of energy has been chosen such that (e) =0, 
with (e") — j deD{e)e^ , and (e^) is a measure of the 
squared width of the band. For the tight-binding band 
ek = 2t{coskx + cos ky) one finds (e^) = 4t^. 

We now use the sum rule on the spectral weight in the 
lower Hubbard band, 



riL 



—ImGiociuj -tO)^hl + S) (B3) 

TT Z 



which, using the analyticity of Gcij-^ — iO) in the lower 
half-plane, is equivalent to the statement 

HL = lim^^oo ujGioc{uj - iO) = lim^^oo I]k 1^ 

= lim^_oof , (B4) 

where is the number of k-points in the first Brillouin 
zone, which are summed over. Combining (B2)-(B4) one 
gets 



1 



(B5) 



The coupling constant / may be related to xioci^) and 
M{uj). Using (11), one may express Gh as 



l'G^,iu;)=M{u;)-xro'M. 



(B6) 



Since Gh is a boson Green's function (see (13)) with posi- 
tive energy spectrum, ujq > 0, the following relation holds 



/•oo 7 

/ —lmGh{uj-tO)^l. 
Jo I" 



(B7) 



From (|B7ll one then finds 
du 



/2 = 



-Im 



M{Lo~iQ)~xJoli^-M)\. (B8) 
A further relation is obtained by using the f-sum rules: 

/cIllj 1 
— ujlmGhi^jJ—iO) — lim uj'^Ghi'^) ^ TT' / 2w„ = cDq 

and 



(B9) 



— XL(^) = lim uj\ioc{oo). (BIO) 

TT cj — >oo 



Taking the limit lim^^^oo ] of ljB6p one finds 

t^q = ^ / —^XlocK^) (Bll) 



with (e^) as defined after (IB2I 
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